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We apply the method of gauge transformation to spin glasses under the microcanonical 
ensemble to study the possibility of ensemble inequivalence in systems with long-range inter- 
' actions and quenched disorder. It is proved that all the results derived under the canonical 

CN ■ 

ensemble on the Nishimori line (NL) can be reproduced by the microcanonical ensemble 
■ irrespective of the range of interactions. This establishes that ensemble inequivalence should 

take place away from the NL if it happens in spin glasses. It is also proved on the NL that 
the microcanonical configurational average of the energy as a function of temperature is ex- 
actly equal to the average energy in the canonical ensemble for any finite-size systems with 
^ ' Gaussian disorder. In this sense, ensembles are equivalent even for finite systems. 
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Equivalence of canonical and microcanonical ensembles is well established for systems 
with short-range interactions.^'^) It is not necessarily the case in the presence of long-range 
interactions because of the absence of additivity:^) If two independent systems with long-range 
interactions are put together, the total energy is not equal to the sum of the energies of two 



^ ■ separate systems, and thus the standard derivation of the canonical ensemble from the mi- 

Q^ , crocanonical ensemble breaks down. One of the prominent consequences is the emergence of 

fN| ' negative specific heat in the microcanonical ensemble, which has long been discussed in the 

astrophysical context*"^) and has been observed in condensed matter of small size where the 
range of interactions is comparable to the system size.^^ Ensemble inequivalence has also been 
discussed in spin systems with long-range interactions with®) and without^"^^) quenched disor- 

•rH , 

^ ■ der. It is therefore important to establish exact /rigorous results on ensemble equivalence and 

H ; 

5t , inequivalence that are applicable generically to a class of problems. The present contribution 

represents a step toward this goal. 

More specifically, we develop a general theory of spin glasses using gauge symmetry in 
the microcanonical ensemble, which should be compared with the corresponding theoretical 
framework in the conventional canonical ensemble, in particular on the so-called Nishimori 
line (NL).i3.i4) ^j^g results agree with those from the canonical ensemble, whenever compar- 
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ison is possible, irrespective of the range of interactions including the infinite-range case. It 
is therefore concluded that canonical and microcanonical ensembles are equivalent in spin 
glasses on the NL as long as physical quantities that can be analyzed by gauge symmetry are 
concerned. 

Let us study the Edwards- Anderson model of spin glasses, 

H = -'^JijSiSj, (1) 

{ij) 

where the Jij are quenched random interactions and the Si{= ±1) are Ising spins. The range 
of interactions is arbitrary and includes the infinite-range Sherrington-Kirkpatrick model. 
The lattice structure and spatial dimensionality are also arbitrary. The present theory is thus 
very general. 

We first discuss the Gaussian distribution of Jij implemented in the microcanonical con- 
text. Our microcanonical analysis starts from the entropy for a given energy, averaged over 
disorder, 

S{E) = ^ I DJ5(iVBJo-5^J.i)logJ]5(i?-i/), (2) 

(ij) s 

where Cj is the normalization factor of the configurational average 

Cj= f BJ6{NBJo-Y.Jij) (3) 

with Nb being the number of bonds, and D J is the Gaussian integral kernel 

In eq. ([2]), the configurational average is taken under the condition Jij — ^bJq- This is 
the microcanonical version of the Gaussian distribution with average Jq, which corresponds 
to the standard (canonical) weight, 

( - ^ Y.^'^^j - -^0?) = ( - ^ + "^0 )) exp (jo J2 -^ij) ■ 

(ij) {ij) (ij) 

The final factor in the above equation indicates that the value of the sum ^ Jij can fiuctu- 
ate according to the probability weight exp(Jo Yl Jij)- The microcanonical version of eq. ([2]) 
strictly enforces the value of ^ Jij to Nb Jo- 
lt is convenient to evaluate here the normalization factor Cj using the condition E = 
—NbJq that will be explained later. 



Cj{E) = j dXexp{-iXE) j DJ expi-iX'^Ji^ 
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= l-jdX exp{-iXE - NBXy2) = ^=i= exp ( - E^ilNB)) . (6) 

Let us apply the gauge transformation Jij — JijUiUj^ Si — Siai where the = ±1 

are gauge variables. Following the prescription established in the canonical ensemble/^' we 
sum the result over all configurations of the gauge variables to have 

SiE) = Wc'jj DJ j;<5(iVBJo-^J^,,c7,a,)log^5(i? + ^J,,.S,5,)), (7) 

cr (ij) S (ij) 

where is the number of sites. The basic strategy of microcanonical calculations is to take 
the derivative of the entropy with respect to the energy and equate the result to the inverse 
temperature f3. If we carefully look at the expression of this derivative, 



^ = aE = Wc]r-' -^J,...^ j:sS(e.e,,j,s,s,)) • 



we notice that the summation over a in the configurational average cancels with the summa- 
tion over S in the denominator when E = —NbJq- This is the microcanonical condition for 
the NL. Then the above expression reduces to 

dS_ _ 

'dE ~ 2^Ci 



13 



^9^5] I BJ5{E + J2J^,S.S,) 
-^fElBJ5{E + Y^J,,) 



(ij) 

dECj{E) E 



(9) 



Cj{E) Nb' 

where we have used eq. ([6]). The second line has been derived by the gauge transformation 
tJ ij y J ij SiS j . Thus we have 

E = -NbI3 (10) 

in precise agreement with the canonical case.^^' The NL condition in the canonical ensemble 
j3 = Jq also follows from the condition E = —Nb-Jq in conjunction with eq. (|10p . It is 
remarkable that both ensembles give the same exact energy as a function of temperature for 
any finite-size system, any dimension, and any range of interactions. Usually, the results of 
two ensembles agree only in the thermodynamic limit. 

Upper and lower bounds for the specific heat C for ii^ = —NbJq can also be estimated as 
shown below, and the results agree with the corresponding canonical inequalities, 

< C < ^. (11) 

Since the specific heat is non-negative and upper-bounded by the same expression as in the 



3/0 



J. Phys. Soc. Jpn. Letter 

canonical ensemble, there exists no ensemble inequivalence under the NL condition at least 
as long as the energy and bounds on the specific heat are concerned. 

To prove eq. (jlip . we evaluate the second derivative of the entropy with respect to the 
energy, using the notation il{E) = J2s ~ ^)-> 



Cj{E) V Cj{E) ) Nb- ^ ^ 



If we remember the relation 

l_ _ d£_ 

it follows 



(13) 



which is the desired inequality pT]) . 

The same analysis applies to the it J model, where Jij = 1 with probability p and Jij = — 1 
with probability 1 — p. We have set J(= \Jij\) = 1 for simplicity without losing generality. 
The points to be modified from the Gaussian case are (i) to replace the integral over Jjj with 
the summation over Jij = ±1, (ii) to replace the normalization Cj with Cp, where 

Cp=Yl ^i^B{2p - 1) - E -^^i) ' (15) 

and (iii) to replace the condition E = —NbJq with E = —Nb{2p — 1). The resulting relation 
corresponding to eq. ([9|) is 

where 

Cp{E) = j (^1^ expi^E + Nb log 2 cosh n). (17) 

Here we have omitted the factor l/27r since it plays no role in the following. The integral 
variable /i corresponds to iX in eq. ([U]). For finite Nb, it is impossible to evaluate the integral 
explicitly, which is a difference from the Gaussian case. In the thermodynamic limit, the 
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saddle-point method yields (3 = fiQ from eq. (fT6]) . where /io is the saddle point specified by 

E + Nb tanh /xq = 0. (18) 

In combination with (3 = fiQ and E = —Nb{2p — 1), we conclude E = —Nstanhf^ under 
the condition tanh/3 = 2p — 1. This is in perfect agreement with the canonical analysis. 
Complete ensemble equivalence holds only in the thermodynamic limit in the it J model. 

Upper and lower bounds for the specific heat on the NL can also be estimated as in the 
Gaussian model. The central inequality is 

dp dlCpjE) ( dEC,{E) V 

dE - Cp{E) V Cp{E) ) ■ ^ ' 

If we naively apply the saddle-point method and take only the leading term, the right-hand 
side reduces fiQ — fiQ = 0, leading to 

1 



p^,<0 (20) 

or C > 0. This positivity of the specific heat is non-trivial in the microcanonical ensemble but 
is not very exciting. A better inequality is obtained from the leading correction to the saddle 
point. The normalization factor Cp{E) is written as, to the leading correction to the saddle 
point, with the notation /(/i) = fJ-E + Nb log 2 cosh /i, 

CpiE) = exp (/(/io)) I d/i exp ((^ - fio?/2a^) , (21) 

where 1/cr^ is the second derivative 9^/(/x)/5/x^|^g = A^'^sech^yUo. The integral converges 
because it runs from —too to too through the saddle point fiQ on the real axis. Similarly, we 
have 

dlCp{E) = exp (/(/io)) j dfxfi^ exp ((/i - /xo)V2ct'). (22) 
Then the right-hand side of eq. p9p is 

dfx exp ((/i - /io)^/2cr^) 

f,l = {-a^ + isl)-f,l = -a\ (23) 



j dfi exp ((/i - fj.of/2a'^] 



It is therefore concluded that 

" CT2 - "iV^sechVo' ^^^^ 
which leads to the bounds for the specific heat < CT^ < A^^sech^/? as is already known in 
the canonical ensemble. 

Identities and inequalities for correlation functions can also be established. Let us realize 
that the microcanonical delta constraint 6{E — H) with E being the control parameter plays a 
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very similar role as the canonical Boltzmann factor exp{—l3H) with /3 the control parameter. 
It is then straightforward to apply the same argument as in the canonical case to derive 
identities and inequalities for correlation functions. If we take the example of the Gaussian 
disorder and use the notation Ej = —NbJq (the NL conditions is E = Ej), the results are 

[((5,5, M"] = K5,5,)s,((5,5,)ijr] (n = l,3,5,---) (25) 

[P{m)] = [Piq]] {E = Ej) (26) 

mS,)E]\ < mS,)Ej\] (27) 

[sgn{SiSj)E] < [sgn{SiSj)Ej], (28) 

where the suffix of the angular brackets specifies the value of the microcanonical energy. The 
square brackets denote the configurational average. The first relation ()25p with n = 1 shows 
that the ferromagnetic correlation on the left-hand side is equal to the spin glass correlation 
on the right-hand side if the system is on the NL. The limit \i — j\ — t- cxd yields m = q, 
where m and q are the ferromagnetic and spin glass (Edwards- Anderson) order parameters, 
respectively. The second identity ()26p restates this fact from the perspective of the distribution 
functions of the magnetization and the spin glass order parameter, 

p. . ^ E5(i),5(^) mQ - si'^S?^)6{E - hW)5{E - g(^)) 

i:si^^,s^^)5{E-Hi^))5{E-Hi^)) ' ^'^^ 

where H^^'^ is the Hamiltonian with the spins Si in eq. ([1]) replaced by the spins 5f ^ of the 
k{= l,2)th replica. This identity (|26p proves that there is no replica symmetry breaking on 
the NL in the sense of non-trivial distribution [P(g)] because the distribution of magnetization 
[P(m)] is always trivial. The third relation ()27p proves that the phase boundary between the 
ferromagnetic and non-ferromagnetic phases below the multicritical point should be either 
vertical or reentrant in the phase diagram. The final inequality ()28p implies that the number 
of mutually parallel spin pairs takes its maximum value on the NL when we change the energy 
(and consequently the temperature). The system is thus in its most ordered state on the NL if 
we focus ourselves to the spin orientation, ignoring the magnitude, as indicated by the signum 
function. All these results are shared by the canonical ensemble. 

The distribution of a local energy —J12S1S2 can be calculated similarly on the NL using 
its gauge invariance. The strategy is exactly the same as in the derivation of eq. Q. The 
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result is 

Me + JuSiS2))ej] = c exp (^-^^^^^(^ _ j^f) , (31) 

where c is the normahzation constant. The local energy naturally distributes in a Gaussian 
form around its mean Jq. 

A different viewpoint can be introduced if we apply random fields and impose a micro- 
canonical constraint that the 'staggered magnetization' along the random fields has a specific 
value. For the Gaussian distribution of disorder, the entropy is 

SiE) = -1- jDJBh6{NBJo-Y,J^j)K^^-Y.^i) 

■ log J2S{E + Y, J^jS^Sj)5[Nm - ^ hiSj) . (32) 

The normalization is now 

Cj= j V>JDh6{NBJo-Y,Jii)^{^h-Y,hi), (33) 

which can be evaluated as in eq. ([6]) under the generalized NL condition E = —NbJq and 
h = m to give 

CjiE,m) = ^^Jj^ exp ( - i?V(2A^i?) " Nm'/2). (34) 

Then, when E = —NbJq and h = m, the first and second derivatives of S with respect to m 
are evaluated as before, 
dS 

— = dmlog C J {E, m) = -mN, (35) 
d'S _ dlCj{E,m) fdmCj{E,m)Y ^^^^ 



dw? Cj{E,m) V Cj{E,m) 

Thus, for m = h > 0, the entropy is a decreasing concave function of the staggered magne- 
tization on the NL (E = —NbJo)- In particular, for m = h = 0, the vanishing value of the 
staggered magnetization is thermodynamically stable in the sense that the entropy is maxi- 
mum. This is natural because the system is not in the spin glass phase (g > 0, m = 0) on the 
NL (q = m) and thus cannot be staggered-magnetized along a given unbiased random field 
(unbiased in the sense /i = 0). For a finite value of h{> 0), the stable value of the staggered 
magnetization is closer to zero than m = h since the entropy will be larger for smaller m 
according to eq. ([55]) . 

In summary, we have shown that all the results obtained from gauge symmetry in spin 
glasses in the canonical ensemble can be reproduced in the micro canonical formulation of 
the same problem. In particular, the micro canonical configurational average of the energy 
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for the Gaussian distribution agrees exactly with the corresponding canonical energy for any 
finite-size systems. This implies complete ensemble equivalence for finite-size systems, an 
unusual phenomenon. These results are valid for a generic system with arbitrary range of 
interactions in arbitrary dimension including the infinite-range limit as long as the system 
is on the NL. We have proved that there are no anomalies such as multiple values of the 
temperature for a given energy or negative specific heat as observed in certain systems with 
long-range interactions.^"^' ^"^^^ If ensemble inequivalence exists in spin glasses with long-range 
interactions, it should happen away from the NL. Preliminary calculations indeed suggest 
possible ensemble inequivalence away from the NL, and the results will be reported in a 
forthcoming publication. 
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